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Abstract 

We consider the stochastic quantization scheme for a non-perturbative stabilization 
of 2D quantum gravity and prove that it does not satisfy the KdV flow equations. It 
therefore differs from a recently suggested matrix model which allows real solutions 
to the KdV equations. The behaviour of the Fermi energy, the free energy and 
macroscopic loops in the stochastic quantization scheme are elucidated. 
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1. Introduction 



There are now at least two proposals for a non-perturbatively stable formulation 
of two dimensional quantum gravity. On the one hand there is the "5th time" 
method [1] rediscovered in the context of Parisi-Sourlas reduction of a one dimen- 
sional supersymmetric matrix model[2] a.k.a. stochastic quantization 1 [4] [5] [6] [7] 

and on the other hand a definition[8] [9] following from non-perturbative KdV 
flows[10][ll]. Both proposals agree[8][2] with the perturbative (i.e. genus by genus) 
results from hermitian matrix models[12] as indeed they must, while also providing 
(apparently) a physically sensible non-perturbative extension. 

By "physically sensible" we mean at least that the string susceptibility and 
correlations of all local operators are real and free of singularities (on the real 
axis). This has been shown to be true in the KdV flow method[8][9], and in the 
stochastic quantization method can be expected from the construction; Recall that 
it is not true of the full non-perturbative result derived from hermitian matrix 
models [13] [14]. We also require that macroscopic loop expectation values be 
physically sensible. For the KdV flow method this has been checked in ref.[9]. 
Part of the purpose of this paper will be to demonstrate similarly sensible results 
for the stochastic quantization method. 2 

Naively one might expect that 2D quantum gravity has a unique non- 
perturbative extension, and hence that the two proposals are equivalent. This 
is not the case. We demonstrate this below, both numerically and analytically. It 
follows unfortunately that at present the question of the "correct" non-perturbative 
extension becomes a philosophical one. The problem is of course that no general 
principle of non-perturbative gravity has yet been defined, and since both methods 
have the same perturbative expansion, it is hard to promote one method over 
another if they have no obvious flaws. 

The proof that the KdV flows are violated non-perturbatively in the stochastic 
quantization scheme will be as follows. The KdV flows induced by the operators 

1 Other suggestions based on stochastic quantization are to be found in ref.[3]. 

2 Our interest in this problem was triggered by some remarks in ref.[15] and it 
seems that our conclusions differ from the the ones presented there. 
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Ok are given by[ll]: 

Here the t^s are couplings conjugate to the operators Ok and p is the string 
susceptibility related to the free energy T by r" = —p. V is given in terms of 
the partition function as Z = exp(r). Primes refer to differentiation with respect to 
p, the cosmological constant. (The string coupling has been absorbed into p and 
the couplings The R^s are the Gelfand-Dikii differential polynomials[16] 

If we assume that no further dimensionful parameter arises at the non- 
perturbative level, then these equations (1.1) are satisfied non-perturbatively if and 
only if p satisfies the string eqn. [9] : 

= pK 2 - ^TZn" + - {K'f (1.2) 

where for the full massive theory 

oo 

n=Y J {k + \)t k Rk[p\ -p . (i.3) 

k=0 

For the purposes of this paper however, we will only be interested in pure gravity 
in which case, without loss of generality, we take only t 2 ^ and 

K = p 2 -p"/Z-n . (1.4) 

Now for negative cosmological constant p, if p satisfies eqn. (1.2) and has an 
asymptotic expansion then it follows from (1.4) and (1.2) that the spherical 
contribution 3 must vanish po = or be imaginary po = ii^ 1 / 2 . (These are the roots 
of the algebraic equation obtained by deleting all derivative terms in (1.2)). For 
negative p an asymptotic expansion for p does indeed exist in the supersymmetric 
ID matrix model since this expansion just corresponds to the WKB expansion. 
Because it is real the leading term is of the form p = K\f\p~\ for some real constant 
k, by scaling arguments. We will prove that k is different from zero. It follows that 
the string susceptibility calculated in the stochastic quantization method cannot 
satisfy (1.2) and hence violates eqns. (1.1) non-perturbatively. 



3 up to analytic terms in p which are always discarded 
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We note in passing that it is in fact natural to introduce a new non-perturbative 
parameter a into (1.2) corresponding to the position of the eigenvalue space 
boundary [17]. It corresponds physically to a world-sheet boundary cosmological 
constant. In this case (1.2) takes the form = (p - a)TZ 2 - \HH" + \ (ll') 2 . It 
is equivalent to (1.2) under the finite 'gauge' transformation exp{— crL_i} and so 
does not alter our conclusions. (eL_i generates the KdV galilean transformation: 
St k = e(k + 3/2)t fc+ i, 5p = 5a = e). 

It may seem a touch paradoxical that non-perturbative information can be 
deduced from perturbative arguments in the regime p < 0. We would therefore like 
to emphasise that the asymptotic results for p < in both schemes should properly 
be regarded as non-perturbative. From the physical point of view both schemes 
derive their connection to pure 2D quantum gravity /0D string theory indirectly 
through exact agreement with the p > hermitian matrix model perturbative 
results, the latter having a direct representation in terms of regularised surfaces 
through the dual of Feynman diagrams. Neither scheme has a direct world- 
sheet interpretation. (This is discussed further for the complex matrix model 
representation of the KdV flow method in ref.[8]). When p < both schemes lose 
their interpretation in terms of surfaces since, if for no other reason, they now 
disagree with the hermitian matrix model. Actually the hermitian matrix model 
has no interpretation in terms of surfaces here either because the p < region is 
outside the radius of convergence of its large iV Feynman diagram expansion. Indeed 
mathematically the perturbative results for p < are separated from those for p > 
by non-smooth behaviour at p = 0. The spherical free energy suffers a 3rd order 
phase transition; In the stochastic quantization scheme for p > supersymmetry 
is preserved to all orders of perturbation theory, while for p < supersymmetry is 
broken even on the sphere. In the representation of the KdV flow scheme in terms of 
an eigenvalue spectrum A G R + the classical eigenvalue distribution hits the "wall" 
at A = as p -> 0+ [8] [9] (more generally[17] A = a as p -> (cr 2 ) + ). 

The rest of this paper is organised as follows: Section 2 derives the spherical 
string susceptibility for p < by calculating the one point function. The role of the 
Fermi energy and the free energy in the super symmetric matrix model is clarified. 
Macroscopic loops on the sphere are discussed in section 3 and a non-polynomial 
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OD hermitian matrix model potential that produces the same results is derived. A 
numerical comparison of the one-point function in both schemes is made in sect. 4, 
and in sect. 5 we draw our conclusions. 

2. The Sphere for Negative Cosmological Constant. 

Recall that the stochastic quantization scheme for the simplest matrix models 
with cubic potentials results in replacing expectation values defined in the hermitian 
matrix model 

< 0(4) >=\jd4 0(cj>) eM-N Tr V(cj>)) (2.1) 

where <p is an N x N hermitian matrix, with expectation values of the same operator 
in the ground state of N non-interacting fermions which have the single particle 
Fokker-Planck (FP) hamiltonian 

Hfp = N[~^ + V fp (X)] (2.2) 

where 

Vfp{\) = \{V'{\)Y- 1 -V"{\) (2.3) 

Whereas the simplest potentials V that give ordinary k = 2 critical behaviour are 
unbounded from below, the Fokker-Planck potential Vpp is not, and this is the 
reason the method leads to non-perturbative stabilization. In fact the expectation 
value of any observable, which in the hermitian matrix model will be given by (2.1) 
, is here defined by 

<0(4>)> F p= j 'dcj> ¥8(0)0(0), (2.4) 

where W (</>) is the groundstate wavefunction of the iV-particle fermionic system. 
In the following we will restrict ourselves to the spherical approximation and to 
observables of the form 0(4>) = Tr f{4>). The spherical approximation is in the 
hamiltonian language equivalent to the WKB approximation and the expectation 
value of the observables Tr /(0) will be given by: 

l<Tr/(0)>=i£rfA jE F -V F p(\)f{\). (2.5) 
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Here A r and A; denote the left- and right turning points, and E F denotes the Fermi 
energy of the system, i.e. the iV'th energy level, scaled by N~ x . The WKB condition 
for the iV'th energy level of the hamiltonian (2.2) is 



1 = I [ d\ JE F - V FP {\) (2.6) 

7T J\ t 



To be explicit let us consider the simplest potential V in the hermitian matrix 
model which leads to k = 2 critical behaviour: 

V(4>) =g<P-<P 3 /3 (2.7) 

The corresponding FP-potential is 

V FP (X;g) = I(£-A 2 ) 2 + A (2.8) 

The critical coupling in the matrix model is g c = 3/2 2 / 3 and the endpoint of the 
eigenvalues distribution is A c = \fg~c/3 = 1/2 x / 3 where V^(A C ; g c ) = g 2 /3 = 3\ 4 and 
the derivatives with respect to A satisfy: V'(X C , : g c ) = V" (A c ; g c ) = 0. 
If we introduce the scaled variables: 

g-g c = g c a 2 ^ \-\ c = \ c ay, a = l/N 2 / 5 (2.9) 

we can write 

V F p(y; n) = V FP (0; y) + X 4 a 3 v fp (y; y) 

3 1 

Vfp{w, v) = -3^y + y 3 + a[ - -/i y 2 + - y A ] 

where V F p(0; fi) is a second order polynomial in g — g c = g c a 2 y. As in [4] the 
FP-hamiltonian can be written 



H FP = N[V FP (0; fi) + \ 4 c a 3 h h 

A 2 

'dy 2 



h f P = - 4 ^- v fp(y^) 



and if we introduce the scaled Fermi energy e/ by 



E F = V F p(0;/j)+X 4 c a 3 e f 
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(2.10) 



the equation (2.6) for the Fermi energy can be written as 

N = 2^ r dy \l e f- v fp(y>ri ( 2 - n ) 



While the original hermitian matrix model even at the spherical level is defined 
only for g > g c , i.e. for \x > 0, the FP-potential is defined for all real \i. Since Ep 
in the spherical limit can be a function of the coupling constant g only (i.e. g c and 
g — g c ), we conclude that 

e/ = |^| 3/2 / (±) (av^). (2-12) 

where f^(-) might be different functions for \x > and /i < 0. Let us assume that 
the /'s have the following expansions 

f (± \x) = ri ±) + 0(x), (2.13) 

where the constants /q are to be determined from (2.11). By differentiating (2.11) 
with respect to \fi\ and taking the limit a — > one gets the following equations 
for/f>: 



0=lim / ay 7 2 . (2.14) 

' /yi V/o Vl 3/2 - "/*(*;/*) 



The integration limits, the turning points, before taking a — > are given by 

y r = V^2: r (/ (±) ) + 0(a), y, = - - - -aa + 0(a 2 ). (2.15) 

Here \J~\/ji\z r denotes the right turning point solution in the limit a — > 0. We have 

/ (±) ±3^-^ = 0. (2.16) 

Subtracting from ±(2.14) the same equation with \i = 0, and then changing variables 
as y — > \f\fji\z we obtain the convergent integrals: 

2Jo - ' dz \-= = --^=] ± Idz 



J-oo . /jt(±) i o„ „3 J- 



V/ (±) ±3z-z 3 V/o (±) ±3z-z 3 



(2.17) 



The solution of (2.17) is well known for fx > due to the inherent supersymmetry 
of the FP-hamiltonian [2], which is unbroken for fi>0 : and is given by 

io (+) = -2 (and z r = -2) . (2.18) 

For fx < it is readily seen that /q - ^ > and z r > 0, but the integrals can not be 
expressed as elementary integrals. Numerical integration gives 

/ (_) « 8.0 (and z r « 1.52) . (2.19) 

Let us now turn to the observables of the theory. The simplest and most 
fundamental ones in the hermitian matrix model are the puncture operator and the 
susceptibility. In the case of the potential (2.7) the 'bare' puncture operator and 
the susceptibitity can be defined as follows 

<°0>=-^5r'( g )=<il\^>. (2.20) 

In (2.20) T denotes the free energy of the hermitian matrix model. The FP-version 
of (2.20) can be found from (2.4) , or in the spherical approximation from (2.5) 

< Oo >fp= - [ X d\ A JE F -V FP {\) = A c (l + [ V dy y Je f -v fp (y) ). 

7T JA, 27T J yi v 

(2.22) 

and we get 

A c a 3 / 2 fvr, y(e' f + 3(y + ±ay 2 )) 



Pfp(v) = --j- d y i — (2-23) 

9c 47T J yi Je f -v fp (y) 



Inserting the integration limits (2.15) and subtracting, as in (2.17), the expression 
for fx = 0, leads to the following result : 

Pfp(v) ~ Pfp(0) = ^ (cie'^ja + 0(( V&) 3 / 2 )) (2.24) 

where 



If we recall that g — g c = g c a 2 fx, the leading non-analytic power of g — g c is given by 
the term involving the Fermi energy. Since it is different from zero for both positive 
and negative cosmological constants, according to (2.18) and (2.19) we conclude that 
the (1.2) must be violated even at the spherical level by the stochastic regularization. 
If we integrate (2.24) with respect to g — g c we get 

< O >fp= y (c + c 2/ ua 2 - e f a s + 0((y^a) 7 / 2 )) (2.25) 

and from (2.24) and (2.25) it is clear that the leading non-analyticity in \x (or g — g c ) 
is entirely contained in the Fermi energy e/ which therefore can be considered as 
the puncture operator 4 . 

Let us end this section with a discussion of the role of relations like (2.20): 
< Oq >= —T'(g)/N 2 in the context of stochastic quantization. Due to the hidden 
supersymmetry, which is unbroken at the spherical level for \i > it is clear that such 
a relation can not be true if we for T use the energy of the iV-Fermi system, since 
the total energy of the iV-Fermi system is exactly zero for all \x > 0. However, it is 
possible to have a somewhat similar formula in the stochastic scheme. By definition 
we have from the second equality in (2.20) that < Oq >fp=< Tr <p/N >fp and 
we can get this quantity from the free energy Tpp of the N- Fermi system if we 
add a source term J Tr (j)/N to the FP-hamiltonian. This term explicitly breaks 
supersymmetry and Tpp(J) will be different from zero and we have, as noticed in 
[6]: 

<0 > FP = dTFP d f 9) \j=o (2.26) 

The modified single particle hamiltonian will be 

d 2 

hfp = -^—+v fp {y,^j,a) (2.27) 
where j is an appropriately scaled source term and 

v fp (y, ^ j, a) = -(Sfi + j)y + y 3 +[ - -^ay 3 + ^y 4 ] (2.28) 

4 This identification of the Fermi energy as the puncture operator is only valid 
for the simplest potentials V(4>) which allow a factorization of the ground state \&o 
in single particle wave functions [18]. 
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It is now clear that taking the limit a — > before calculating expectation values 
means that we are breaking supersymmetry in a way identical to adding a source 
term j since only the combination 3fj, + j appears in (2.28) after a is zero. Of course 
the integrals defining the expectation values will now be divergent since we have 
dropped the stabilizing term in the potential. However, if one just introduces a 
cut-off A ~ —1/a in the integrals it can be shown that the leading non-analytical 
behaviour is still the same as with the full potential. In this way one can convert 
(2.26) to an equation similar to (2.20) where the differentiation is with respect to 
\x and where Tpp(p,) is the energy calculated from the iV-Fermion system with the 
unbounded potential — 3[iy + y 3 , but with a lower cut-off A ~ — 1/a. All results in 
this section can be derived using this technique. 

3. Macroscopic Loops on the Sphere. 

Motivated by ref.[15], we make some remarks on the expectation value of a 
single macroscopic loop in the stochastic quantization scheme. We only consider 
the spherical limit, however as we have seen above the behaviour of the spherical 
limit for \i < is in a sense non-perturbative. It is not unreasonable to expect the 
full non-perturbative result to have a similar qualitative (but more smoothed out) 
behaviour. 

From identification through the (dual of the) Feynman diagram expansion in 
the hermitian matrix model we know that the operator representing the insertion 
of an n- polygonal loop is given by 



This class of observables belongs to the ones which can readily be defined in the 
stochastic regularization scheme. In the spherical approximation (2.5) we have 



For finite n each observable clearly has the same critical behaviour as the puncture 
operator. However interesting behaviour arises if one scales n — > oo as 



<W(n) >=< ^- Tt(/) n > 



(3.1) 




(3.2) 



n = l/a. 
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(3.3) 



In this way we can identify / with a macroscopic length of the loop, a being identified 
with the "lattice" spacing. Introducing the scaled variables (2.9) in (3.2) we get 



< W{ a ] >= if 27 ly dV ^f^VfM (1 + *y) l/a - (3-4) 
In the scaling limit (1 + ay) l l a — > e ly and we get 

/ \ l / a 1 rvi i i 

< W(~) >= -jj- < >; < w(l) >= — Jjy je f -v fp (y)e l y. (3.5) 

This expression is perfectly well defined for both positive and negative \i as long as / 
is positive. For / negative (3.1)-(3.4) are not well defined and this is reflected in the 
exponential divergence of (3.5) for a — > 0. As is well known [13] [19] [20] the Laplace 
transform of w(l) is the generator of one point functions and has very transparent 
analyticity properties. It is defined by 



F W = M e" k < u,(l) > = J- f% ^'' V,M (3.6) 



o 2tt J yi y-z 

From this integral representation it is clear that F(z) is analytic everywhere in the 
complex plane except along the cut z G [yi, y r ] and we can express < w(l) > by the 
inverse laplace transform: 

— e zl F(z) . (3.7) 

-ioo+r Zni 

Here r is to the right of any singularities of F(z) in the z complex plane, i.e. to the 
right of y r . In fact (3.7) makes sense even for negative / and can be considered as a 
definition of negative length loops. For / negative we can clearly close the contour 
to the right in (3.7) and we get 

< w(l) >= for / < . (3.8) 

This physically very reasonable result is true independently of whether ji is positive 
or negative 5 . 



5 We should note that our definition seems to differ from the one given in [15] 
where one gets an expectation value of w(l) for / < 0. 
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Up to now the there has been no difference between negative and positive 
cosmological constants \i. However the cut-structures of 2nuf p (z) = \Jef — Vf p (z) 
in the complex plane are different. The hidden super symmetry, unbroken for /i > 0, 
ensure that in this case e/ — Vf p (z) has a double zero yo = \f\ii\+ 0(a) and in this 
case (recall that yi = —A/a + 0(a)): 

2iru fp (z) = (y - z)yj (y r - z)(z - yi)a/A . (3.9) 

However, for [i < the double zero yo splits in two complex zeroes given approxi- 
mately by 



Z± =-^-± t ^-^y2 + 4\^\ (3.10) 

and 

2nu fp (z) = \J(z- z+)(z - z-)(y r - z)(z - yi)a/A (3.11) 

has an additional cut from z+ to Z-. From the integral representation (3.6) it is 
clear that F(z) falls off as 1/z for \z\ — > oo, is analytic except along the cut [yi, y r ] 
on the real axis and that 

F(y + ie)-F(y-ie) =2niu(y). (3.12) 

When u(z) is given by (3.9) the cut of u(z) coincide with the one of F(z) and 
F(z) is uniquely determined as 

F(z) = ^V'(z) + i7ru(z) (3.13) 

where ^V'(z) is given by (2.7) and determined uniquely by the requirement that 
F(z) = 0(l/z) for \z\ ->oo. 

For \i < the additional cut in u(z) invalidates (3.13), and one has to subtract 
a function with this cut and which falls of as 1/z. It is uniquely determined as 

i X '(z) = ^-i dz'^1 (3.14) 
2 K ' 2ixi Jc z' - z v ' 

where the contour C surrounds the new cut. F is now given as 

F(z) = \(V'(z) + X'(z)) + mu(z). (3.15) 
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From the classic work [21] of Brezin et al. it is known that if the eigenvalue 
distribution of the OD-matrix model in the large N limit is given by u(y), defined 
on an interval [j//, y r ], then the resolvent F(z), defined by 6 

F(z) = / dy , 

J vi y~ z 

is related to the potential V e ff(<f>) of the matrix model by 

F(y±ie) = ^V' eff (y)±iiru(y), ye[y t ,y r }. (3.16) 

From (3.13) and (3.15) we conclude that the OD matrix model effective potential 
which gives the same results as the FP-potential in the spherical limit is given by 

V eff = V + e(-n)X . (3.17) 

Note that X is non-polynomial. It would be interesting to examine whether the full 
stochastic quantization scheme solution would also follow from V e //. 



4. A Numerical Comparison. 

We have shown analytically in sects. 1 and 2 that the stochastic quantization scheme 
is different from the KdV flow method, by comparing their leading asymptotics as 
H — > — oo. It is of interest to compare them in the small [i regime too. We can do 
this numerically using the one-point function for the puncture operator, < Oo >, 
calculated in ref. [6] in the stochastic quantization scheme and <Oq> calculated in 
the KdV flow method in ref. [9]. 

We refer the reader to these papers for more details on those calculations. In 
ref. [6] the real and imaginary parts of the non-perturbative hermitian matrix model 
solution are also computed and plotted. We have copied the five points for the 
stochastic quantization scheme (the circles) and for the real part of the hermitian 

6 We have the relation 

1 1 a 3/2 
R( w )= <Tr- >= F(z), w = X c (l + az) 

JM <p — W A c 

where R(w) denotes the resolvent in unsealed variables. 
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matrix model solution (the triangles) onto fig.l, using the conventions of ref.[6]. 
The remaining line in fig.l is calculated in the KdV flow method [9]. The errors 
in this line are negligable (~ 10~ 5 ). Although the three graphs clearly differ, it is 
intriguing that they are qualitatively so similar. 

The translation from the conventions of ref. [9] to those of ref. [6] is as follows: 

P-Po = 8v / 3^ 3 / 2 (<Oo>-^ 3/2 ) 

is the sphere subtracted one-point function, up to factors as indicated, and is called 
T and Q in ref. [6] . The horizontal coordinate z is proportional to the cosmological 
constant: z = \(2\f?>) 2 ^ p. These translations are derived using the recursion 
relations of the orthogonal polynomial coefficients f n and a n (eqns.(26) and (27) of 
ref. [6]) to calculate the spherical contribution of ref.[6]'s definition of the one point 
function 

1 1 N ~ 1 

< O0> .<_Tr^>=— £a„ 

in the contour rotated hermitian matrix model. Using the continuum limits given 
there in eqns.(29-31) 7 we find 



(I- 9/ 9c) 3 / 2 

where / satisfies the Painleve I equation in the form 

d 2 f 



dy f(y) (4.1) 



48z 5 / 2 dy 2 



+ f z = y 



We use this and eqn.(1.4) to make the fundamental identifications / = —p{x)/^JJi, 
y = x/fi and z = |(2y / 3) 2 / 5 ^t from which it follows that (4.1) translates to 



8^ 



dx p(x) 



7 Eqn.(30) should read a 2 y n = ■ ■ ■ 
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5. Conclusions 

As shown above the KdV flow symmetry of 2D quantum gravity is broken 
non-perturbatively by the real solution provided by stochastic quantization scheme. 
Since another scheme exists which preserves these flows non-perturbatively and 
which allows for a real solution [9], and since we have found no obvious unphysical 
behaviour of physical observables in the two schemes (as discussed in sec. 3), one 
is left with ambiguities. On the one hand we have the general method of stochastic 
stabilization of bottomless actions [1] , which works in any dimension. As already 
pointed out in ref. [6] it disagrees with the stabilization via contour rotation [22] 

. In 2d gravity the contour rotation method satisfy the KdV-flow symmetry, 
but leads to unacceptable complex solutions. On the other hand the unique real 
solution to (1.2) is an alternative candidate for a non-perturbative definition of 2d 
gravity [9]. The preference of this definition to the one of stochastic stabilization 
would amount to the declaration that KdV-flows are a fundamental property of non- 
perturbative 2D quantum gravity. Keeping in mind on the one hand the ultimate 
goals of understanding higher dimensional gravity and string theory, and on the 
other, the lack of understanding of the fundamental physics behind the KdV-flow 
structure, it seems to us rather that the results obtained in this article highlight the 
need for a better understanding of the basic principles underlying non-perturbative 
quantum gravity. 
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Figure Captions 



Fig. 1: The expectation value of the puncture operator with the sphere contribu- 
tion removed; z is proportional to the cosmological constant. The triangles 
are the real part of the hermitian matrix model result. The circles are 
values given by the stochastic quantization scheme, and the unmarked 
line is the result of the KdV flow scheme. 



17 



